Shape resonances for ultracold atom gases in carbon nanotube waveguides 
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We propose an experimentally viable setup for the realization of one-dimensional ultracold atom 
gases in a nanoscale magnetic waveguide formed by two doubly-clamped suspended carbon nano- 
tubes. All common decoherence and atom loss mechanisms are shown to be small. We discuss 
general consequences of a non-parabolic confinement potential, in particular novel two-body shape 
resonances, which could be observed in this trap. 

PACS numbers: 03.65.Nk, 03.75.Gg, 73.63.Fg 



The ongoing progress in the fabrication and ma- 
nipulation of micro- or nanoscale structures has re- 
cently allowed for systematic studies of ultracold atom 
gases, where current-carrying wires generate magnetic 
fields trapping neutral atoms ('atom chips') 0, @- 
For instance, the Bose-Einstein condensation (BEC) of 
microchip-confined atoms has been successfully demon- 
strated by several groups ||. Such an approach is par- 
ticularly promising in the context of integrated atomic 
matter- wave interferometry and optics Q , and combines 
the strengths of nanotechnology and atomic physics. So 
far, decoherence and atom loss constitute central impe- 
diments, since atoms are relatively close to 'hot' macro- 
scopic surfaces or current-carrying wires (with typical 
diameters of several fJJXi), where the Casimir-Polder po- 
tential and Johnson noise can seriously affect stability 
[EL IE 0- To reduce these effects, further miniaturiza- 
tion would be desirable. While at first sight this goal 
conflicts with the requirement of large currents forming 
the trapping potentials, we propose that when using sus- 
pended carbon nanotubes (NTs) Q (with diameters of 
a few nm) as wires, nanoscale atom chip devices basi- 
cally free of decoherence or atom loss can be built with 
state-of-the-art technology. With relevant length scales 
below optical and cold-atom de Broglie wavelengths, this 
also paves the way for the observation of interesting and 
largely unexplored many-body physics in one dimension 
(ID) |9j, involving either bosons or fcrmions. Examples 
include the interference properties of interacting matter 
waves [TJJ . spin-charge separation 

[HI 

, fractional statis- 
tics and atom number fractionalization [12| . and the ID 
analogue of the BEC-BCS crossover . Previous real- 
izations of ID cold atoms were reported using optical lat- 
tices 0,0,0] and magnetic traps 0, but they involve 
arrays of ID or elongated 3D systems, where the above 
many-body effects are difficult to observe. Our proposal 
completely eliminates unwanted substrate effects and im- 
plies a drastically reduced transverse size (a few nm) of 
the cloud. Hence rather high atom densities could be 
achieved, allowing for stable operation and sensible opti- 
cal detection schemes. 

The proposed nanoscale waveguide confining ultracold 
atoms to ID is sketched in Fig. ^ The setup employs 



two suspended doubly-clamped NTs, where nanofabrica- 
tion techniques routinely allow for trenches with typical 
depth and length of several /jm Q. To minimize deco- 
herence and loss effects [fj, the substrate should be in- 
sulating apart from thin metal strips to electrically con- 
tact the NTs. Since strong currents in the mA-regime 
are necessary, thick multiwall nanotubes (MWNTs) or 
'ropes' (bundles) Q are best suited. Due to the sus- 
pended geometry, the disturbing influence of the sub- 
strate is largely eliminated, and employing an additional 
longitudinal magnetic field to eliminate Majorana spin 
flips El, 

neutral atoms in a weak-field seeking state 
can be trapped between the NTs. Studying various 
sources for decoherence, heating or atom loss, and es- 
timating the related time scales, we find that, for rea- 
sonable parameters, detrimental effects are small. We 
also analyze effects due to the non-parabolicity of the 
transverse trap potential, which is inevitable in all com- 
mon traps but particularly pronounced in the present 
case. We find many two-body shape resonances, which 
allow to tune atom-atom interactions similar to standard 
Feshbach resonances. Previous work has only studied 
parabolic traps, where center-of-mass (COM) variables 
decouple. Then incoming scattering waves can visit only 
one out of many bound states normally available, leading 
to a single 'confinement-induced resonance' [2(1 • Here we 
solve the two-body problem for general transverse con- 
finement, and then apply the results to the proposed trap. 
As a concrete example, we shall consider 87 Rb atoms in 
the weak-field seeking hypcrfine state \F,mp) — |2, 2). 

We next describe the setup in Fig. ^ in detail, where 
the same (homogeneous) current I flows through the two 
NTs at (±#0,0, z), and a magnetic field B z is applied 
along the z-direction. Neglecting boundary effects due 
to the finite tube length L, the magnetic field at position 
x = (x,y,z) = (x ± ,z) is 
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with the vacuum permeability fj,Q. The transverse con- 
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finement potential is V(xjJ = /x|B(x)|, where \i = 
mFgFfJ-B with the Lande factor gp and the Bohr mag- 
neton [Xb- Under the adiabatic approximation mp 
is a constant of motion, and the potential is harmonic 
very close to the minimum of the trap, with frequency 
lu = [fj,/ '{mBz)] 1 / 2 [i I / '(ttXq) and associated transverse 
confinement length a±_ = {h/muj) 1 / 2 <C xq, where to is 
the atom mass. The adiabatic approximation is valid as 
long as w < Wi with the Larmor frequency ujl — pB z /h. 
Trapped atoms can in general also make non-adiabatic 
Majorana spin flips to a strong-field seeking state, and 
thereby escape from the trap The associated loss 

rate is T ~ (ttw/2) exp(l - 1/x) with X = huj/{fiB z ) flU ]. 
For convenience, we switch to a dimensionless form of the 
full potential V(xj_) by measuring energies in units of hu> 
and lengths in units of ax, 
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which depends only on d — xo/a± and x- To give an 
example, we take 1 = 1 mA, representing a reason- 
able current through thick MWNTs @. For d = 10 
and T/oj = 10~ 6 (\ = 0.067), we obtain B z = 20 G, 
xq = 79 nm and uj = 2tt x 1.85 MHz, corresponding to a 
very tight trap. The potential is shown for these param- 
eters in Fig. [3 HU . 

For stable operation, it is essential that destructive ef- 
fects like atom loss, heating or decoherence are small. 
One loss process proceeds via (i) non-adiabatic Majorana 
spin flips as discussed above. Another one is (ii) atom 
loss due to tunneling out of the trap. The WKB tun- 
neling rate -ft can be estimated easily for an atom in the 
transverse ground state escaping along the least-confined 
y-direction, sec Fig. [5J 



7t/w ~ (2tt) 1 exp 
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For the above parameters, we find numerically y\ = 
1.47,2/2 = 68.09, and hence j t /uj w 10~ 100 . Atom loss 
may also originate from (iii) noise-induced spin flips, 
where current fluctuations cause a fluctuating magnetic 
field generating the Majorana spin flip rate 
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At room temperature and for typical voltages Vq ~ 1 V, 
we have Hujl <C ksT <c eVo, and Si {ujl) is expected 
to equal the shot noise 2e//3 of a diffusive wire. For 
the parameters above, a rather small escape rate re- 
sults, 7 s f « 0.017 Hz. Next we study (iv) the trans- 
verse NT deflection due to their mutual magnetic re- 
pulsion, using a standard elasticity model for a dou- 
bly clamped wire in the limit of small deflections [23j |. 
For small NT displacements <p{z), the equation of mo- 
tion is pl4> = —YMi(f)"" + p I 2 /{Attxq), where pi is 



the linear mass density, Y is Young's modulus and Mj 
the NT's moment of inertia. The static solution un- 
der the boundary conditions 0(0, L) = (()'{0,L) = is 
4>{z) = hq[Iz{z — L)] 2 / {96nY MiXq) . Using typical ma- 
terial parameters from Ref. @, the maximum displace- 
ment is <f>{L/2) 0.08nm for L = 10/im. Hence the 
mutual magnetic repulsion of the NTs is weak, (v) Ther- 
mal NT vibrations might create decoherence and heat- 
ing, and could even cause a transition to the first excited 
state of the trap. The maximum mean square displace- 
ment is a 2 = ((f> 2 {L/2)} = k B TL 3 /(l92YMi) 23], which 
for the above parameters gives a w 0.2 nm at room tem- 
perature. This is much smaller than the transverse size 
a± of the atomic cloud. Detailed analysis shows that 
the related decoherence rate is also negligible. Another 
decoherence mechanism comes from (vi) current fluctua- 
tions in the NTs. Following the analysis of Ref. 0, the 
corresponding decoherence rate can be estimated for the 
above parameters as 7 c /w < 10~ 7 . We conclude that 
no serious decoherence, heating or loss mechanisms are 
expected for reasonable parameters of this nanotrap. 

Next we discuss s-wave atom-atom interactions in this 
trap. We assume that the 3D interaction can be de- 
scribed by a Fermi pseudopotential [24( via the 3D scat- 
tering length a. To keep generality, let us investigate two- 
body scattering for arbitrary confining potential V(xj_), 
and later specialize to Eq. (JTJ. Unlike for a harmonic 
trap, the problem does not decouple when formulated in 
terms of the COM coordinate R = (xi+x 2 )/2 = (Rj_, Z) 
and the relative coordinate r = x 2 — xi = (rj_, z). Follow- 
ing standard steps [jj , the pseudopotential is enforced as 
a boundary condition for the two-particle wavefunction, 



*(R,r 



o) = t?t(i-M/<0- 
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Without loss of generality, we set the longitudinal COM 
momentum to zero such that Z drops out. Furthermore, 
we can put rj_ = and then let z — > in Eq. @ ■ The so- 
lution of the two-particle Schrodinger equation ^(Rj_, z) 
then takes the general form 



*(R±, z) = *o(R±, z)+ JdR' ± G E (B.±, z; R x , 0) 



TO 

(3) 

For ^ = 0, this leads to bound-state solutions with 
energy E < 2E discussed elsewhere j2^, where Eq 
is the single-particle ground-state energy. Here we fo- 
cus on scattering solutions at low energies E slightly 
above 2Eq, where exactly one transverse channel is open. 



Then * (Rx,2) 



Akz „/,2 



^o(R-l) describes two incom- 



ing atoms with (small) relative longitudinal momentum 
hk = y/2m{E — 2Eq) in the (transverse) single-particle 
state V'o(R-l) with energy Eq. Some algebra yields the 
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two-particle Green's function 



G E (R±, z; R±, 0) = ^(R±)^(R!l) 



im 
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where the bar denotes complex conjugation, and 
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Since / e (fc) = —1/(1 + tom) follows from Eq. JSJl, 
aiD can be identified with the ID scattering length. 
The effective ID atom-atom interaction potential is then 
Vm(z, z') = girv5(z — z') with interaction strength giD = 
— 2fi 2 /(maiD) lU H3- For very low energies, k — ► 0, 
we can now put £ = 2£?o in Eq. (|1U|) . For a binding 
trap, C,2E is an Hcrmitian operator with discrete spec- 
trum {A„} and eigenvectors |e„), and hence we find 



hD - 2c 2 E 



l(0|e, 



A n + l/(47ra) 



(11) 



Here, G\ ^ is the (transverse) single-particle Green's func- 
tion, with single-particle states tp\ an d energy E\, and 
the two-particle Green's function Gt acts on the Hilbert 
space 7i c iosod, orthogonal to the open channel. Inserting 
Eq. J3J into Eq. J3J yields for \z\ — > oo a standard scat- 
tering solution *(R, z) = V> 2 (RjJ(e ifcz + / e (/fc)e ifc l z l) with 
scattering amplitude 



f e (k) = - / dR^ 2 (Rl)/(R4 



(5) 



Enforcing Eq. (J2J then leads to an integral equation, 

1 = J dR' x ( E (R ± , R' ± )/(R1) (6) 



47ra 



-V>o 2 (Rj-) 
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dR!4 2 (Rl)/(Rl), 



where the singular behavior has been split off by intro- 
ducing a regularized kernel in H c i OS ed) 



Ce = 



dt 



(e £t G t (R ± ,R' ± ) - £^(R± - R±) 



The integral equation © is most conveniently solved by 
expansion in a suitable orthonormal basis \j), 

l/> =X)/iU>, /j = / dR±(j|R±)/(R±), (7) 



where |0) corresponds to (RjjO) = c^q(Rj_) with nor- 
malization constant c. Thereby we can express Eq. © 
in compact notation, 



Ana c 
which is solved by 

-l/c 
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Here, the parameter aiD follows in the form 

2c 2 
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(0|[C B + l/(4™)]- 1 |0)- 
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Let us then denote -ffciosed as the projection of the Hamil- 
tonian to W c k>sed- The boundary condition for the corre- 
sponding bound state is enforced by — | /) / (47ra) = ( E \f) ■ 
Thus, for a = — (47rA„) , a bound state of -ffciosed 
with energy E — 2E exists, fulfilling Eq. (0) with 
f(R) = (Rj_|e n ). In agreement with general arguments 
|2(ll2fij. provided (0|e„) ^ 0, a bound state of -ffciosed at 
a = — (47rA„) _1 corresponds to a shape resonance in the 
ID interaction strength giD. 

The two-body problem can thereby be solved numer- 
ically for any given confinement potential by diagonal- 
izing (,2E - We note in passing that for certain poten- 
tials, including a hard-box confinement, Eq. Ill 111 can 
also be computed analytically in closed form |2fil |. In 
general, every eigenvalue A„ corresponds to a different 
shape resonance, unless the overlap (0|e„) vanishes due 
to some underlying symmetry. For a parabolic con- 
finement, the decoupling of the COM motion implies 
that only one resonance is permitted, and some alge- 
bra gives indeed from Eq. Hll|) the known result aiD = 
— (a 2 L /a)(l — 1.0326a/aj_) [2fJ. For a general confine- 
ment, there are in principle infinitely many resonances. 
In practice, however, only few of them can be resolved 
since most of them appear at a/a± — > 0. In addition, 
most resonances are extremely sharp and will therefore 
be hard to detect. The ID interaction strength #id for 
the potential (JIJ in the case of 87 Rb atoms is shown in 
Fig. |21 as a function of the NT current I. At least two 
different shape resonances are predicted to be observable 
for realistic and discernible values of I in the mA regime. 
These resonances could be used to tune the strength and 
the sign of the two-body interactions by simply adjusting 
the NT current. 

To conclude, we propose a nanoscale waveguide for ul- 
tracold atoms based on doubly clamped suspended nan- 
otubes. All common sources of imperfection can be made 
sufficiently small to enable stable operation of the setup. 
Detection certainly constitutes an experimental challenge 
in this truly ID limit. However, we note that single- 
atom detection schemes are currently being developed, 
which would also allow to probe the tight ID cloud here, 
e.g., by combining cavity quantum electrodynamics with 
chip technology |2j, or by using additional perpendicu- 
lar wires/tubes 'partitioning' the atom cloud |22j. This 
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may then allow to study interesting many-body physics 
in ID in an unprecedented manner. Atom-atom interac- 
tions can be tuned by shape resonances, which we have 
described for arbitrary transverse confinement. 
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FIG. 2: (Color online) Transverse trapping potential of the 
nanoscale waveguide, see text. 
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FIG. 3: (Color online) The solid curve gives gm as function 
of I for the potential JJ| via numerical solution of Eq. pip 
for d ~ 13 and \ — 0.067. The corresponding parabolic pre- 
diction for the same u is shown as dashed curve. 



